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CHAPTER 1 



The CFL phase and m s : An effective field theory approach 



Thomas Schafer 



Department of Physics 
North Carolina State University 
> ! Raleigh, NC 27695 

We study the phase diagram of dense quark matter with an emphasis on the role 
of the strange quark mass. Our approach is based on two effective field theories 
(EFTs). The first is an EFT that describes quark quasi-particles near the Fermi 
surface. This EFT is valid at energies small compared to the chemical potential. 
The second is an EFT for the Goldstone modes in the paired phase. We find that 
in response to a non-zero strange quark mass the CFL phase first undergoes a 
transition to a kaon condensed phase, and then to a gapless phase with a non-zero 
Goldstone boson current. 



o 

1. Introduction 



Searching for exotic states of matter at high baryon density and high temperature 
is one of the central efforts in nuclear and particle physics. Calculations based on 
weak-coupling QCD indicate that the ground state of thre e flav or baryonic matter 



at very high density is the color-flavor- locked (CFL) phase HEED The CFL phase is 
characterized by a pair condensate 

WC-ystf) = - SfStW. (1) 
This condensate leads to a gap in the excitation spectrum of all fermions and com- 
pletely screens the gluonic interaction. Both the chiral SU(3)l x SU(3)r and color 
SU(3) symmetry are broken, but a vector-like SU(3) flavor symmetry remains un- 
broken. 

At baryon densities relevant to compa ct st ars distortions of the ideal CFL state 
due to quark masses cannot be neglected GEO The most important effect of a non- 
zero strange quark mass is that the light and strange quark Fermi momenta will 
no longer be equal. When the mismatch is much smaller than the gap one finds for 
degenerate quarks, we expect that it has little consequence, since at this level the 
original particle and hole states near the Fermi surface are mixed up anyway. On 
the other hand, when the mismatch is much larger than the gap, we expect that 
the ordering one finds in the symmetric system is disrupted, and that to a first 
approximation one can treat the light and heavy quark dynamics separately. 



1 
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This argument is qualitatively right, but the correct picture turns out to be 
much more complicated, and much more interesting. The phase diagram of cold 
dense strange quark matter contains phases with kaon and eta condensates, phases 
with non-zero currents, and crystalline states. Performing systematic calculations 
of the properties of these states is not straightforward, even if the coupling is weak. 
A standard set of tools are Dyson-Schwinger equations. This approach becomes 
quite involved once the strange quark mass is taken into account, because there 
are many more gap p arame ters, and maintaining electric neutrality and color gauge 
invariance is difficult 11 ' V However, since chiral symmetry is broken in the CFL 
phase we know that the dependence on the quark masses is constrained by chiral 
symmetry. It is therefore natural to study the problem using effective field theories. 
In practice we will employ a two-step procedure. In the first step we match the 
microscopic theory, QCD, to an effective field theory of quasi-particles and holes 
in the vicinity of the Fermi surface. In the second step we match this theory to an 
effective chiral theory for the CFL phase. 

2. High density effective theory 

2.1. Effective field theory near the Fermi surface 

The QCD Lagrangian in the presence of a chemical potential is given by 

c = $W + wo-m)i>-\g%,g%,, (2) 

where — + igA^ is the covariant derivative, M is the mass matrix and fi is 
the baryon chemical potential. If the baryon density is very large perturbative QCD 
calculations can be simplified. The main observation is that the relevant degrees of 
freedom are particle and hole excitations in the vicinity of the Fermi surface. We 
shall describe these excitations in terms of the field ib v (x), where v is the Fermi 
velocity. At tree level, the quark field -0 can be decomposed as -0 = Vv+ + VV- 
where ib v ,± = Pv,±4> with P„,± = i(l ± a. ■ v)ip. Note that P v .± is a projector on 
states with positive/negative energy. To leading order in we can eliminate the 
field ib- using its equation of motion. The lagrangian for the ibx. field is given by 

i^rrnrm 

C = ^^iv.D-^- 9 -^^y v -lG^G^ + .... (3) 

with w M = (l,v). Note that v labels patches on the Fermi surface, and that the 
number of these patches grows as u 2 . The leading order v ■ D interaction does not 
connect quarks with different v, but soft gluons can be exchanged between quarks in 
different patches. In addition to that, there are four, six, . . . fermion operators that 
contain fermion fields with different velocity labels. These operators are constrained 
by the condition that the sum of the velocities has to be zero. 

In the case of four-fermion operators there are two kinds of interactions that 
satisfy this constraint, see Fig. ^ The first possibility is that both the incoming 
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Fig. 1. Kinematics of four-fermion operators in the effective theory. 

and outgoing fermion momenta are back-to-back. This corresponds to the BCS 
interaction 

c = 4 E v f r ' R f r '^ • (V:,rvi„,), (4) 

v',r,r' 

where v-v 1 — cos 6 is the scattering angle, Rf r (x) is a set of orthogonal polynomials, 
and r, r' determine the color, flavor and spin structure. The second possibility is 
that the final momenta are equal to the initial momenta up to a rotation around 
the axis defined by the sum of the incoming momenta. The relevant four-fermion 
operator is 

C = \ £ FF'RF' (v • (Vlr'4). (5) 

In a system with short range interactions only the quantities f}(0) are known as 
Fermi liquid parameters. 

The effective field theory expansion is complicated by the fact that the number 
of patches N v ~ /i 2 /A 2 grows with the ratio of chemical potential n over the EFT 
cutoff A. This implies that some higher order contributions that are suppressed by 
l//i 2 can be enhanced by powers of N v . The natural solution to this problem is to 
sum the leading order diagrams in the large N v limit ^1 In the gluon sector this 
corresponds to summing particle-hole loops in gluon n-point functions. There is a 
simple generating functional for thes e loop integrals which is known as the hard 
dense loop (HDL) effective actional 

m 2 ^ v a v^ b 

Lhdl- — G ^ (v-DY ^' ^ ' 

In perturbation theory the dynamical gluon mass is given by m 2 = Nfg 2 [i 2 /(4ir 2 ). 
In the fermion sector we have to consider four-fermion operators at leading order. 
Operators with six or more fermion fields are suppressed. The effective lagrangian 
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is given by 

£ = v4 Uv-d- ^ - \g;„g; u + c HDL (7) 



U 2 S l 



r L 



2.2. Non- Fermi liquid effects 

In this Section we briefly study the effective field theory in the regime uj < m 
where u> is the excitation energy and m is the effective gluon mass H^l. Since electric 
fields are screened the interaction at low energies is dominated by the exchange of 
magnetic gluons. The transverse gauge boson propagator is 

A 3 (fc)= '^T^L ' (8) 
M - k 2 + i^m?^S- 

U A k 

where we have assumed that |A:o| < We observe that the propagator becomes 
large in the regime \ko\ ~ \k\ 3 /m 2 . If the energy is small, |fco| m, then the typical 
energy is much smaller than the typical momentum, 

|fc| ~ (m^fcol) 1 ^ » |fc |. ( 9 ) 

This implies that the gluon is very far off its energy shell and not a propagating 
state. We can compute loop diagrams containing quarks and transverse gluons by 
picking up the pole in the quark propagator, and then integrating over the cut in 
the gluon propagator using the kinematics dictated by equ. 10 . In order for a quark 
to absorb the large momentum carried by the gluons and stay close to the Fermi 
surface this momentum has to be transverse to the momentum of the quark. This 
means that the term k 2 _/(2n) in the quark propagator is relevant and has to be kept 
at leading order. Equation Q shows that fc^/(2/i) ^> ko as ko — > 0. This means 
that the pole of the quark propagator is governed by the condition fc|| ~ k\/(2a). 
We find 

fc^«? 2/ V /3 fco /3 , fc||~ 5 4 V /3 fco 2/3 . (io) 

In the low energy regime propagators and vertices can be simplified even further. 
The quark and gluon propagators are 

= ', — 7' (") 

Po~P\\- 2^ +iesgn(p ) 

Dij(k)= , ~ lSij i , (12) 

and the quark gluon vertex is gvi(\ a /2). Higher order corrections can be found by 
expanding the quark and gluon propagators as well as the HDL vertices in powers 
of the small parameter e = (ko/m). 
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We will refer to the regime in which all momenta, including external ones, satisfy 
the scaling relation QlOjl as the Landau damping regime. The Landau damping 
regime is completely perturbative, i.e. graphs with extra loops are always suppressed 
by extra powers of e 1 / 3 . Every quark propagator scales as e -2 / 3 , gluon propagators 
scale as e -2 / 3 , and every loop integral gives e 7 ^ 3 . The quark-gluon vertex scales as 
e° and the HDL three-gluon vertex scales as e 1 / 3 . Using these results we can show 
that additional loops always increase the power of e 1 / 3 associated with the diagram. 

The effective theory describes a non-Fermi liquid. This is clear from the ap- 
pearance of fractional powers and logarithms in the low energy expansion. The 
simplest diagram which gives a logarithmic term is the fcrmion self energy. We find 

nafanns 

The scale inside the logarithm is determined by matching the effective theory in 
the Landau damping regime to an effective theory that contains electric gluon ex- 
changes. The po log(po) term leads to a vanishing Fermi velocity as we approach the 
Fermi surface. 

The effective theory can also be used to study other higher order corrections. 
We find, in particular, a QCD version of Migdal's theorem: In the Landau damping 
regime loop corrections to the quark-gluon vertex are suppressed by powers of e 1 / 3 . 



2.3. Color superconductivity 

It is well known that the particle-particle scattering amplitude in the BCS chan- 
nel q(p) + q{—p) — > q(p') + q{—p') is special. The total momentum of the pair 
vanishes and as a consequence loop corrections to the scattering amplitude are log- 
arithmically divergent. This implies that all ladder diagrams have to be summed. 
Crossed ladders, vertex corrections, etc. are perturbative and follow the scaling rules 
discussed in the previous section. 

If the interaction in the particle-particle channel is attractive then the BCS 
singularity leads to the formation of a pair condensate and to a gap in the fermion 
spectrum. The gap can be computed by solving a Dyson-Schwinger equation for 
the anomalous (particle-particle) self energy. In QCD the interaction is attractive 
in the color anti-triplet channel. The structure of the gap is simplest in the case of 
two flavors. In that case, there is a unique color anti-symmetric spin zero gap term 
of the form 

WC 75 ^) ~ Ae 3ab e«. (14) 
Here, a, b labels color and i, j flavor. The gap equation is given 

by |19|2U|21|22| 

A(p0 =*[ d qi lo g ( ^ - A(g4) , (15) 
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where the scale Abcs = 2567r 4 (2/iV/) 5 / 2 <7~ 5 /i is again determined by electric gluon 
exchanges. The solution to the equation was found by Son El The value of the gap 
on the Fermi surface Ao = A(po = 0) is 

A ~2A BCS exp(-^^^exp(-^j . (16) 

This result is correct up to 0(g) corrections to the pre-exponent. In order to achieve 
this accuracy the q 2 ui Iok(u>) term in the normal self energy, equ. {TSJl, has to be 
included in the gap equation | 23 | 24 | 25 | rp^g condensation energy is given by 

where N4 = 4 is the number of condensed species. 

The situation is slightly more complicated in QCD with Nf = 3 massless flavors. 
One possibility is to embed the Nf = 2 order parameter into Nf = 3 QCD. This 
option is usually called the 2SC phase. The energetically preferred phase is the CFL 
phase described by the order parameter given in equ. (JJJ. In the CFL phase there 
are eight fermions with gap Acfl and one fermion with gap 2Acfl- The value of 
the gap in the CFL phase is smaller than the one in the 2SC phase by a factoJ2 
2-1/3 rpjjg con( j ensa y orl energy in the two phases is given by 

6 = (i^) { 12 • 2- 2 / 3 CFL (18) 
and the CFL phase is preferred by a factor (27/4) 1 / 3 ~ 1.89. 



2.4. Mass terms 

Mass terms modify the parameters in the effective lagrangian. These parameters 
include the Fermi velocity, the effective chemical potential, the screening mass, the 
BCS terms and the Landau parameters. At tree level the correction to the Fermi 
velocity and the chemical potential are given by 

v F = l-— T , on = -- — . 19) 

The shift in the Fermi velocity also affects the coupling gv F of a magnetic gluon to 
quarks. It is important to note that at leading order in g this the only mass correction 
to the coupling. This is not entirely obvious, as one can imagine a process in which 
the quark emits a gluon, makes a transition to a virtual high energy state, and then 
couples back to a low energy state by a mass insertion. This pro cess would give an 
0(mj fi) correction to <?, but it vanishes in the forward direction ESI. 

Quark masses modify quark-quark scattering amplitudes and the corresponding 
Landau and BCS type four-fcrmion operators. Consider quark-quark scattering in 
the forward direction, v + v' — * v + v'. At tree level in QCD this process receives 
contribution from the direct and exchange graph. In the effective theory the direct 
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Fig. 2. Mass terms in the high density effective theory. The first diagram shows a O(MM^) term 
that arises from integrating out the field in the QCD lagrangian. The second diagram shows 
a 0(M 2 ) four-fermion operator which arises from integrating out tp— and hard gluon exchanges. 



term is reproduced by the collinear interaction while the exchange terms has to be 
matched against a four-fermion operator. The spin-color-flavor symmetric part of 
the exchange amplitude is given by 

M (v,v';v,v') = -I^-^ll-^T^-l (2°) 
v ; AN c N f p F \ p F l-xf y ' 

where Cp — — 1)/(2N C ) and x = v ■ v' is the scattering angle. We observe 
that the amplitude is independent of x in the limit m — > 0. Mass corrections are 
singular as x — > 1. The means that the Landau coefficients Fi contain logarithms 
of the cutoff. We note that there is one linear combination of Landau coefficients, 
Fo — i*i/3, which is cutoff independent. 

Equations H19I20J) are valid for Nf > 1 degenerate flavors. Spin and color anti- 
symmetric BCS amplitudes require at least two different flavors. Consider BCS 
scattering v + (—v) — > v' + (-V 1 ) in the helicity flip channel L + L — > R + R. The 
color-anti-triplet amplitude is given by 



C A g 2 m 1 m 2 
4 V F Pf 



M(v,-v;v,-v) = —-s-—s—. (21) 



where m\ and mi are the masses of the two quarks and Ca — {N c + l)/(2N c ). We 
observe that the scattering amplitude is independent of the scattering angle. This 
means that at leading order in g and m only the s-wave potential Vq is non-zero. 

In order to match Green functions in the high density effective theory to an 
effective chiral theory of the CFL phase we need to generalize our results to a 
complex mass matrix of the form C = —^lMiPr — TpRM'ipL, see Fig. [5] The 5/j, 
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term is 



— (i>l + MM^ L+ + ^ R+ M^Mi> R+ 



2p F 



(22) 



and the four-fermion operator in the BCS channel is 

-,2 



c 



9 



(^C^mC^){Tr [X A M(X D fX B M(X c f] (23) 
- ilr [X A M(X D ) T ] Tr [A B M(A C ) T ] }. 



Here, we have introduced the CFL eigenstates ip A defined by ip® = i/j A (X A ) a i/V^, 
A = 0,...,8. 



2.5. Normal quark matter 

Before we consider quark mass effects in the superconducting phase we would like to 
review mass corrections in the normal phase of quark matter. We will consider three 
flavor quark matter with massless up and down quarks and massive strange quarks. 
There are some simple mass effects that we can directly deduce from equ. (|19(l . The 
strange quark chemical potential is shifted by 5fj, s = m 2 /(2pp) with respect to the 
strange quark Fermi momentum, and the strange quark Fermi velocity is smaller 
than one. 

It would seem that perturbative corrections to these results are sensitive to mo- 
menta far away from the Fermi surface and cannot be computed in the framework 
of an effective field theory for modes near pp. Landau showed, however, that Galilei 
invariance leads to a relation between the chemical potential and the Fermi mo- 
mentum that only involves the interaction on the Fer mi surface. This relation was 
generalized to relativistic systems by Baym and Chin 1221. They showed that 



pdji = 



PF 



dpt 



(24) 



where Fq 1 are color-flavor-spin symmetric Landau coefficients. We observe that 
this relation involves precisely the combination of Landau coefficients that does 
not depend on the EFT cutoff. This implies that we can integrate equ. (I24|l as 
in ordinary Landau Fermi Liquid theory. The re sult can be used to determine to 
thermodynamic potential to 0{g 2 ). We getGHEH 



12tt 2 



1) 



—to H — to m 
2 ' 2 



16tt 3 



2, M 

i In — 



TO 

2 



2u q 



(25) 



with u = \J /-J.' 2 — m? . This result corresponds to a momentum space space subtrac- 
tion scheme. The thermodynamic potential in the MS scheme was computed by 
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9 



Fraga and Romatschke Perturbative corrections reduce the pressure of a quark 
gas. The same is true for mass corrections to the pressure in the non-interacting 
system. The 0(a s m 2 ) term, however, increases the pressure. This is related to the 
fact that the exchange energy in a degenerate quark gas changes sign in going from 
the non-relativistic limit, dominated by the Coulomb interaction, to the relativistic 
system in which magnetic interactions are more important. 

In applications to neutron stars we have to include weak interactions and enforce 
electric charge neutrality. Weak interactions can convert strange quarks into up 
quarks, electrons and neutrinos. We will assume that neutrinos can leave the system. 
Under these assumptions the system is characterized by a baryon chemical potential 
/i and an electron chemical potential The electron chemical potential is fixed by 
the condition of electric charge neutrality, (dfl) / (dfi e ) = 0. 

To leading order in rn 2 s /p 2 F the electron chemical potential is given by 



At tree level we find /i e = /i s /2 and the Fermi surfaces are split symmetrically 



Perturbative corrections reduce the magnitude of u e . In fact, since the 0(a s ) term 
is enhanced by a large logarithm log(pF/rn), the electron chemical potential can 
become negative. This result is probably not reliable. In particular, the 0(a s ) term 
is modified if the MS mass is used. 

3. Chiral theory of the CFL phase 
3.1. Introduction 

The main topic of this review is the effect of the strange quark mass on the CFL 
phase. In principle this problem can be studied by solving a gap equations which 
includes mass corrections to the chemical potential, the Fermi velocity, and the 
quark-quark interaction. In full QCD this has not been attempted yet, but there 
are many calculations of this type that are based on effective four-fermion models 
of QCD. We will compare our results with some of these calculations in Sect. 13.81 

There are several difficulties with microscopic calculations of this kind. The 
first is that they typically require an ansatz for the gap parameter. We will see 
that both flavor symmetry and isospin are broken, and that the gaps for the left 
and right handed fermions acquire a relative phase. This means that the number 
of independent gap parameters is quite large. The second difficulty is that while 
the ideal CFL state is automatically neutral with respect to all charges, this is no 
longer the case once the state is perturbed by a finite strange quark mass. This 
means that we have to introduce flavor chemical potentials and gluonic background 
fields. Finally, we know that chiral symmetry is broken in the CFL phase and this 
implies that physical observables are non-analytic in the quark mass. 




(26) 



Pf=Pf + /^A Pf = Pf-M s /2- 



(27) 
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Fig. 3. Contribution of the 0(M 2 ) BCS four-fermion operator to the condensation energy in the 
CFL phase . 



These non-analyticities are related to Goldstone bosons. This suggests that we 
should study the effects of a non-zero strange quark mass using an effective field 
theory for the Goldstone modes. Indeed, it is well known that chiral symmetry places 
important constraints on the mass dependence of QCD observables, and that these 
constraints are most easily implemented by using an effective lagrangian. 



3.2. Chiral effective field theory 

For excitation energies smaller than the gap the only relevant degrees of freedom are 
the Goldstone modes associated with the breaking of chiral symmetry and baryon 
number. Since the pattern of chiral symmetry breaking is identical to the one at 
T = ii = the effective lagrangian has the same structure as chiral perturbation 
theory. The main difference is that Lorentz-invariancc is broken and only rotational 
invarian ce is a good symmetry. The effective lagrangian for the Goldstone modes is 
given by "2 

f 2 

C efJ = ^Tr [VoSVoS 1 - ^<9;E<9 i E t ] + [STr(MSt) + h.c] 
+ [AiTr(ME t )Tr(MS t ) + A 2 Tr(ME t ME + ) 

+A 3 Tr(ME t )Tr(M t E) + h.c] + .... (28) 

Here E = exp(i<j> a X a / /„.) is the chiral field, f n is the pion decay constant and M is a 
complex mass matrix. The chiral field and the mass matrix transform as E — ► LYiR^ 
and M — ► LMR^ under chiral transformations (L, R) € SU(3)l X SU(3)r. We have 
suppressed the singlet fields associated with the breaking of the exact U(l)v and 
approximate U(1)a symmetries. 

At low density the coefficients / T , B, Ai, . . . are non-perturbative quantities that 
have to extracted from experiment or measured on the lattice. At large density, on 
the other hand, the chiral coefficients can be calculated in perturbative QCD. For 
the derivative terms this is most easily done by matching the two-point functions 
of flavor currents between QCD and the chiral theory. The mass terms can be 
computed by matching the mass dependence of the vacuum energy. At leading 
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order in ot s the Goldstone boson decay constant and velocity are^^ 

21-81og(2) f&\ 2 = 1 
U 18 \2n 2 )' 3' 1 ' 

Mass terms are determined by the operators studied in Sect. 12.41 We observe that 
both equ. 11221) and Il2i'>[l are quadratic in M. This implies that B — in perturbative 
QCD. B receives non-perturbative contributions from instantons, but these effects 
are small if the density is large, see Sect. 13.71 

We observe that X L = MM^ /{2p F ) and X R = M^M/(2p F ) in equ. (|22H act as 
effective chemical potentials for left and right-handed fermions, respectively. For- 
mally, the effective lagrangian has an SU(3)l x SU(3)r gauge symmetry under 
which Xl r transform as the temporal components of non-abelian gauge fields. 
We can implement this approximate gauge symme try in the CFL chiral theory by 
promoting time derivatives to covariant derivatives' 



/MAft\ (M^M\ 



(30) 



The BCS four-fermion operator in equ. (|23| cont ributes to to the condensation 
energy in the CFL phase, see Fig. EI We find I ^ I ^O 

A£ = | (l^M]) 2 - Tr A/ 2 | + <-> M f ) . (31) 

This ter m can be matched against the Ai terms in the effective lagrangian. The 
result i s E32El 



3A 2 

A 1 = -A 2 = A 3 = 0. (32) 

The vacuum energy also receives contributions from Landau-type four-fermion op- 
erators, but these terms are proportional to Tr[MAf^] and do not depend on the 
chiral field E. 

We can now summarize the structure of the chiral expansion in the CFL phase. 
The effective lagrangian has the form 




a) { Pf aJ { P % J 

Loop graphs in the effective theory are suppressed by powers of 9/(4^/^). Since 
the pion decay constant scales as f n ~ p F Goldstone boson loops are suppressed 
compared to higher order contact terms. We also note that the quark mass expansion 
is controlled by m 2 /(p F A). This is means that the chiral expansion breaks down if 
m 2 ~ pfA. This is the same scale at which BCS calculations find a transition from 
the CFL phase to a less symmetric state. 
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Fig. 4. Phase structure of CFL matter as a function of the effective chemical potential fi s = 
m,g/(2pp) and the lepton chemical potential [iq, from Kaplan & Reddy (2001). A typical value of 
fi s in a neutron star is 10 MeV. 



3.3. Kaon condensation 

Using the chiral effective iagrangian we can now determine the dependence of the 
order parameter on the quark masses. We will focus on the physically relevant case 
m s > m u = nid- Because the main expansion parameter is m 2 /(pF^) increasing 
the quark mass is roughly equivalent to lowering the density. The effective potential 
for the order parameter is 



f 2 r 2 

V eff = ^Tr [2X L ZX R tf - X 2 - X 2 R ] - A x [(Tr(M£t)) - Tr ((MS*) 



(34) 



The first term contains the effective chemical potential fi s = m 2 /(2pF) and favors 
states with a deficit of strange quarks (with strangeness S = —1). The second 
term favors the neutral ground state X = 1. The lightest excitation with positive 
strangeness is the K° meson. We therefore consider the ansatz E = exp(iaA4) which 
allows the order parameter to rotate in the K° direction. The vacuum energy is 



V(a) 



■i / 1 1 \ 2 

1 / mt — to 



2 Pl 



sm(aY + {m° K y (cos(a) - 1) , 



(35) 



where {m° K ) 2 — {AA\ / f 2 )m{m + m s ). Minimizing the vacuum energy we obtain 



cos(a) 



1 n s < m° K 

(m° K ) 2 . 



The hypercharge density is 



1 «) 4 
A* 4 



(36) 



(37) 



February 9, 2008 0:32 



WSPC/Trim Size: 9.75in x 6.5in for Review Volume 



int'bcs 



The CFL phase and m s 



13 



m 




ITU 



Fig. 5. Phase structure of CFL matter as a function of the light quark mass m and the strange 
quark mass m s , from Kryjevski, Kaplan & Schafer (2005). 

This result has the same structure as the charge density of a weakly interacting 
Bose condensate. Using the perturbative result for A\ we can get an estimate of the 
critical strange quark mass. We find 



m s (crit) =3.03-™y 3 A 2/3 , 



(38) 



from which we obtain m s (crit) ~ 70 MeV for A ~ 50 MeV. This result suggests 
that strange quark matter at densities that can be achieved in neutron stars is 
kaon condensed. We also note that the difference in condensation energy between 
the CFL phase and the kaon condensed state is not necessarily small. For fi s — > A 
we have sin (a) — > 1 and V(a) — > / 2 A 2 /2. Since / 2 is of order h 2 /(2tt 2 ) this is 
comparable to the condensation energy in the CFL phase. 

The strange quark mass breaks the 5(7(3) flavor symmetry to SU(2)i x U(l)y- 
In the kaon condensed phase this symmetry is spontaneously broken to U(l) o. If 
isospin is an exact symmetry there are two exactly massless Goldstone modes ™sl 
the K° and the K + . Isospin breaking leads to a small mass for the K + . The phase 
structure as a function of the strange quar k m ass and non-zero lepton chemical 
potentials was studied by Kaplan and Reddy ESI se e Fig. 0] We observe that if the 
lepton chemical potential is non-zero charged kaon and pion condensates are also 
possible. 



3.4. Eta meson condensation 

The CFL phase also contains a very light S — mode which can potentially become 
unstable. This mode is a linear combination of the rj and rj and its mass is propor- 
tional to m u rrid- Because this mode has zero strangeness it is not affected by the 
fi s term in the effective potential. However, since m u ,rrid <C m s this state is sen- 
sitive to pert urb ative a s m 2 corrections. The relevant contribution to the effective 
lagrangian is 



SC = -8A [(TrMX^) 2 + Tr(M£ t ) 2 ] + h.c. 



(39) 
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Fig. 6. This figure shows the fermion spectrum in the CFL phase. For m B = there are eight 
fermions with gap A and one fermion with gap 2 A (not shown). Without kaon condensation 
gapless fermion modes appear at fi s = A (dashed lines). With kaon condensation gap less modes 
appear at fi s = 4 A/3. 



with 

M-^AJ, A? = «.^. (40, 



Here, Aq is a color-flavor symmetric gap parameter which is generated by perturba- 
tive corrections to the dominant, color-flavor anti-symmetric gap^. Since the one- 
gluon exchange interaction in the color-symmetric channel is repulsive the 0(a s m^) 
contribution to the mass of the 7/ — if mode tends to cancel the 0(m u m,i) term. 
When the two terms become equal the eta condenses. The resulting phase diagram 
is shown in Fig.0 The precise value of the tctra-critical point (m* , m*) depends sen- 
sitively on the value of the coupling constant. At very high density to* is extremely 
small, but at moderate density m* can become as large as 5 MeV, comparable to 
the physical values of the up and down quark mass. 



3.5. Fermions in the CFL phase 

So far we have only studied Goldstone modes in the CFL phase. However, as the 
strange quark mass is increased it is possible that some of the fermion modes become 
light or even gapless In order to study this question we have to include fermions 
in the eff ective field theory. The effective lagrangian for fermions in the CFL phase 

is EEEi 



C = Tr (Nkv^DpN) - DTr (TV V 7g N yj _ FTl (n V 7g ^ N ^ 

+ T { ( Tr ( NlNl ^ ~ [Tr (iVi)]2 ) -i L ^ R ) + h - c ] ■ ( 41 ) 
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Nl.r are left and right handed baryon fields in the adjoint representation of flavor 
SU(3). The baryon fields originate from quark-hadron complementarity I^QI. We 
can think of N as describing a quark which is surrounded by a diquark cloud, 
Nl ~ <1l(qlQl)- The covariant derivative of the nucleon field is given by D^N = 
d^N + i[Vfj,, N]. The vector and axial-vector currents are 

= -\ + ^d^} , A, = - (V^Et) £, (42) 

where £ is defined by £ 2 = S. It follows that £ transforms as £ — ► L£ i U(x) Jf = 
U(x)£ltf with U(x) G SU(3) V . For pure SU{3) flavor transformations L = R = V 
we have U(x) = V. F and D are low energy constants that determine the baryon 
axial coupling. In perturbative QCD we find D = F = 1/2. 

The effective theory given in equ. I|41|l can be derived from QCD in the weak 
coupling limit. However, the structure of the theory is completely determined by 
chiral symmetry, even if the coupling is not weak. In particular, there are no free pa- 
rameters in the baryon coupling to the vector current. Mass terms are also strongly 
constrained by chiral symmetry. The effective chemical potentials (Xl,Xr) appear 
as left and right-handed gauge potentials in the covariant derivative of the nucleon 
field. We have 

D N = d N + i[T ,N], (43) 
To - ~\ {do + iX R ) £t + £+ (d + iX L ) £} , 

where X L = MM^ /(2pp) and X R = M^M/(2p F ) as before. (X L ,X R ) covariant 
derivatives also appears in the axial vector current given in equ. 1)42(1 . 

We can now study how the fermion spectrum depends on the quark mass. In 
the CFL state we have £ = 1. For /i s = the baryon octet has an energy gap A and 
the singlet has gap 2 A. As a function of /i s the excitation energy of the proton and 
neutron is lowered, u> Ptn = A — /x s , while the energy of the cascade states 5 _ ,S° 
particles is raised, cjh = A + fi s . All other excitation energies are independent of 
fj, s . As a consequence we find gapless (p,ri) and (S _ ,S°) _1 excitations at ^ s = A. 

The situation is more complicated when kaon condensation is taken into ac- 
count. In the kaon condensed phase there is mixing in the (p, £ + ,£~,S~) and 
(n, E°, S°, A 8 , A ) sector. For m° K « (i s « A the spectrum is given by 

t+f^' <*n*B?A={ A, (44) 

4Ms ' [ 2A. 

Numerical results for the eigenvalues are shown in Fig. H3 We observe that mixing 
within the charged and neutral baryon sectors leads to level repulsion. There are 
two modes that become light in the CFL window fj, s < 2A. One mode is a linear 
combination of proton and E + particles, as well as S~ and E~ holes, and the other 
mode is a linear combination of the neutral baryons (n, E°, S°, A 8 , A ). 




9, 2008 0:32 



WSPC/Trim Size: 9.75in x 6.5in for Review Volume 



int'bcs 



16 T. Schdfer 



/'- 




Fig. 7. Left panel: Energy density as a function of the current jk for several different values 
of ii s = m^/(2pp) close to the phase transition. Right panel: Ground state energy density as 
a function of fi s . We show the CFL phase, the kaon condensed CFL (KCFL) phase, and the 
supercurrent state (cur KCFL). 



3.6. Meson supercurrent state 

Recently, several groups have shown that gapless fermion modes lead to instabilities 
in the current-current correlation functionE^GS Motivated by these results we have 
examined t he stab ility of the kaon condensed phase against the formation of a non- 
zero current [ 43 | 44 | c ons icler a spatially varying U(l)y rotation of the maximal kaon 
condensate 

/l 

U(x)£, K oU\x) = 1/V2 ie**W/y/2\. (45) 
\(He- ife ( x V\/2 l/\/2 

This state is characterized by non-zero currents 



(46) 



In the following we compute the vacuum energy as a function of the kaon current 
j K = V(])k- The meson part of the effective lagrangian gives a positive contribution 





















ie 


-i<t>K 






£ = \vlfhl- (47) 

A negative contribution can arise from gapless fermions. In order to determine this 
contribution we have to calculate the fermion spectrum in the presence of a non-zero 
current. The relevant part of the effective lagrangian is 

C = Tr (NUv^D^N) + Tr (V75 (p A + v ■ A\ N\ 

+ ^ {Tr (NN) - Tr (N) Tr (N) + h.c.} , (48) 
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where we have used D = F = 1/2. The covariant derivative is D N = doN+i[py, N] 
and DiN = diN + iv ■ [V, N] with V, A given in equ. (gnjl and 

The vector potential py and the vector current V are diagonal in flavor space while 
the axial potential pa and the axial current A lead to mixing. The fermion spectrum 
is quite complicated. The dispersion relation of the lowest mode is approximately 
given by 

. (l-l ) 2 3 1_ _ 
^ + 2A ~4 Ms ~4 Jk ' (50) 
where / = v ■ p — pf and we have expanded loi near its minimum Iq = (p s + v ■ 
jV)/4. Equation l|5U|l shows that there is a gapless mode if p, s > 4A/3 — Jk/3. The 
contribution of the gapless mode to the vacuum energy is 



where dfl is an integral over the Fermi surface. The integral in equ. (|51() receives con- 
tributions from one of the pole caps on the Fermi surface. The result has exactly the 
same structure as the energy functional of a non-relativistic two-component Fermi 
liquid with non-zero polarization, see Sect. 21 Introducing dimensionless variables 

3JL 3p s - 4A 



we can write £ = cAffh{x) with 

f h (x)=x 2 -- \{h + xfl 2 <d{h + x)-{h-xfl 2 Q{h-x) . (53) 
x L 

We have defined the constants 



C : 



«=T^J 5 (54) 



where c w = (21 — 81og(2))/36 is the numerical coefficient that appears in the weak 
coupling result for /„.. According to the analysis inEHlthe function fh(x) develops 
a non-trivial minimum if hi < h < hi with h\ ~ —0.067 and hi ~ 0.502. In 
perturbation theory we find a = 0.43 and the kaon condensed ground state becomes 
unstable for (A - 3/^/4) < 0.007A. 

The energy density as a function of the current and the groundstate energy 
density as a function of p s arc shown in Fig. In thes e pl ots we have included 
the contribution of a baryon current jb, as suggested in EH. In this case we have 
to minimize the energy with respect to two currents. The solution is of the form 
Jb ~ 3k- The figure shows the dependence on jk for the optimum value of jb- We 
have not properly implemented electric charge neutrality. Since the gapless mode 
is charged, enforcing electric neutrality will significantly suppress the magnitude of 
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Fig. 8. Screening mass of flavor gauge fields in the CFL phase. The two curves show the second 
derivative of the effective potential with respect to the current at the origin and at the minimum 
of the potential. 

the current ESI We have also not included the possibility that the neutral mode 
becomes gapless. This will happen at somewhat larger values of fj, s . 

We note that the ground state has no net current. This is clear from the fact 
that the ground state satisfies S£ / 'S(V '4>k) — 0. As a consequence the meson current 
is canceled by an equal but opposite contribution from gapless fermions. We also 
expect that the ground state has no chromomagnetic instabilities. From the effective 
lagrangian we can compute the screening length of a SU(3) p gauge field. In the CFL 
phase the isospin and hypercharge screening masses are 



which shows the magnetic instability for h > and has the characteristic square root 
singularity observed in microscopic calculations In Fig. [8] we show the screening 
mass in the kaon condensed CFL phase and the supercurrent phase as a function 
of fi s . We observe that there is an instability in the homogeneous phase, but the 
instability disappears in the supercurrent state. 

3.7. Instanton effects 

The results discussed in Sects. I?"3l3.6l are based on an effective chiral theory of the 
CFL phase. The coefficients of the effective theory arc determined by matching to 
perturbative QCD calculations. We find that the chiral coefficients are "natural", 
i.e. their numerical value agrees with simple dimensional estimates. For example, 
f% is given, up to a coefficient of order one, by the density of states on the Fermi 
surface. This suggest that many of our results are valid even if the QCD coupling- 
is not weak. 




(55) 
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Fig. 9. Instanton contribution to the BCS four-fermion operator in QCD with three flavors. 



A possible exception is the linear mass term in the effective chiral theory. The 
coefficient B in equ. I|28l) vanishes to all orders in perturbation theory, but B re- 
ceives non-perturbative contributions from instantons. In QCD with three flavors 
instantons induce a four-fermion operator 

(2N C -1 



C 



2(2^)V 



{$R,h° >V'Z,,g 1 )(V'fl,/2 fT ^L,g 2 ) + [M <r+ M^L <r+ R)\ , (56) 



1 

~8N~ C 

see Fig- El Here, fi,gt are flavor indices and a^ v = §[7^,71,]. The instanton size 
distribution n(p, p) is given by 



n{p, p) = C N 

Cn — 



/8tt 2n 


2N C 


r stt 2 1 


U 2 , 


) P~ 5 exp 


i g(p) 2 \ 



exp [—Nfp 2 p 2 



0.466 exp(- 1.679^)1.34^ 



(57) 
(58) 



(AT c -l)!(AT c -2)! 

At zero density, the p integral in equ. l|5fi[) diverges for large p. This is the well-known 
infrared problem of the semi-classical approximation in QCD. At large chemical 
potential, however, large instantons are suppressed and the typical instanton size 
is p ~ p^ 1 -c A-qcd wnere Aqczj is the QCD scale parameter. The instanton 
contribution to the vacuum energy density is 

-1 2 



S = - 



16 



i(p,v)dpY( 7T P) 4 P 3 



9 



A 



2tt 2 



Tr [M + M f ] 



(59) 



This result can be matched against the O(M) term in the effective chiral lagrangian. 
We find 

2 



B = C 



N " 



8tt 4 r(6) 
3 3 6 



3V27T f£_ 



9 



V 2tt 2 



8tt^ " 
9 2 



A 



QCD 
p 



A 



QCD' 



(60) 



where we have performed the integral over the instanton size p using the one- 
loop beta function. The coefficient B is related to the quark-anti-quark conden- 
sate, (rpip) = —2B. The linear mass terms contributes to the mass of the kaon, 
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Table 1. Effective chemical potential for the different quark modes. The quark 
modes are labeled by their color (rgb) and flavor (uds). We also show the Q 
charge and the corresponding baryon mode. The leading order result in the 
neutral CFL phase corresponds to /i e = fig = and fig = — 



mode 


Q 


effective chemical potential 


leading order 


baryon mode 


ru 





-|(Ue + ^M3 + J|U8 


MO 




gd 





+ 3^ e — 5^ 3 + 3^ 8 


MO 


(Ao,A s ,S ) 


bs 





i 1 2 
+ g/i e — 3M8 — Ms 


MO 




rd 


-1 


+ §Me + |/«3 + |M8 


MO 


£- 


gu 


1 


2 1,1 
— 3^ e — 2^ 3 + 3^ 8 


MO 


E+ 


rs 


-1 


+ l^ e + 1^3 + §M8 - Ms 


MO - Ms 




bu 


1 


-§^ e - §M8 


MO + Ms 


P 


gs 







MO - Ms 


H° 


bd 





1 1 2 


MO + Ms 


n 



8m 2 K ~ (4B / f%)mm s , and tends to inhibit kaon condensation. At moderate density 
this contribution is quite uncertain because the result is very sensitive to the value of 
the strong coupling constant. Using the one-loop running coupling and Aqcd — 200 
MeV gives values as large as Suik ~ 100 MeV. This is almost certainly an overesti- 
mate because the main contribution comes from large instantons with size p ~ 0.5 
fm. If the average instanton size is constrained to be less than the phenomenological 
value at zero density, p = 0.35 fm, then we find Suik < 10 MeV. 



3.8. Microscopic models 

In this section we shall compare some of our results with micros copic calculations 
based on Nambu-Jona-Lasinio (NJL) models, see for example l ^7 | 48 | ^ natural 
choice is to model the interaction between quarks by a local four-fermion inter- 
action with the quantum numbers of one-gluon exchange 

£ = 4, (ip + M7o — M) ip + G (V^A» 2 ■ ( 61 ) 

This lagrangian has the symmetries of QCD except that the SU (3) color symmetry 
is global rather than local, and the axial U (1)a is not anomalous. The effects of the 
anomaly can be taken into account by adding the instanton induced four-fermion 
interaction given in equ. (|56|l . Since color is a global symmetry the NJL model does 
not exhibit a Higgs mechanism and color superconductivity leads to an extra octet 
of colored Goldstone bosons. 

The pairing ansatz is usually taken to be 

(CC75V>)> ~ A l£ Qb V + A 2 e a <% 2 + A,/"'"',,,,. (62) 

This ansatz reduces to the CFL ansatz for Ai = A 2 = A3 and allows for the 
breaking of hypercharge and isospin. The ansatz is not sufficiently rich to allow for 
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kaon condensation. The chiral field £ can be denned as 



S = XY\ 

where X, Y are related to the left and right handed condensates 



(63) 



(64) 
(65) 



In order to study kaon condensation we have to allow the left and right handed 
order parameters t o be in dependent. This problem was studied in two papers by 
Buballa and Forbes 

Color and electric charge neutrality are enforced by adding chemical potentials 
(//3, na, n e ) for color isospin, color hypercharge and electric charge. In most micro- 
scopic calculations the effect of the strange quark mass is taken into account by 
considering an effective chemical potential /i s = m s /(2pp) for the strange quark. 
In Table^we list the chemical potentials for the different quark states. The quarks 
are labeled by color (rgb) and flavor (uds). They are grouped into the main pairing 
channels in the CFL phase. We also show the leading order result in the charge neu- 
tral CFL phase. In this case we find /i e = A*3 = and fis = ~l J -s and the dominant 
pair breaking stress occurs in the (rs) — (bu) and (gs) — (bd) sector. 

We can associate the quark states in the microscopic theory with the baryon 
fields in the effective theory by computing their quantum numbers under the un- 
broken SU(3)f symmetry. The effective theory is formulated in terms of gauge 
invariant fields and does not require color chemical potentials. In weak coupling we 
can derive the effective theory by integrating out gluonic degrees of freedom. In this 
case the equations of motion for the color gauge potential automatically enforce 
color neutrality. 

We observe that the leading order results in the NJL calculation agree with the 
results derived from the effective field theory. The comparison was extended to the 
kaon condensed phase by Forbes He finds that the bu (proton) and bd (neutron) 
states split, and that the pair breaking stress on both states is reduced. The critical 
Hs for gapless states is [i s — 1.2 A, in qualitative agreement with the leading order 
EFT result pt s = 4A/3. 

Microscopic calculations show that in the CFL phase the gap parameters Ai,2,3 
are very similar even if the pair breaking stress p, s is close to the value of the gap 
S3. The main effect is a reduction of the average gap which is presumably related 
to the decrease in the common Fermi momentum that occurs if fj,s is increased at 
constant /is. In the gapless CFL phase, on the other hand, the splitting between 
the different gap parameters becomes significant. Neither one of these two effects 
is included in the EFT calculations discussed in Sects. 13.313.61 The response of the 
gap to the trace part of MM^ jpp can be taken into account by adding a higher 
order operator to the effective lagrangian equ. (|4*T|) . but the coefficient of this term 
is model dependent. Indeed, in perturbative QCD the gap will likely increase if pf 
is lowered. The splitting between the gaps in the gapless CFL is due to the fact that 
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Fig. 10. Conjectured phase diagram for a polarized cold atomic Fermi g function of 

the scattering length a and the difference in the chemical potentials Sfj, = ah — AH' from Son & 
Stephanov (2005). 

gapless regions on the Fermi surface no longer contribute to pairing. In the EFT 
treatment this effect is not included because the current is assumed to be much 
smaller than the gap. Clearly, in the regime where the number of gapless states is 
large the back- reaction on the gap has to be taken into account. 

4. Cold atomic systems 

A nice system in which stressed pairing can be studied in the laboratory is a 
cold dilute gas of fermionic atoms. Using Feshbach resonances it is possible to 
tune the interaction between the Bose-Einstein (BEC) limit of tightly bound di- 
atomic molecules and the Bardeen-Cooper-Schrieffer (BCS) limit of weakly corre- 
lated Cooper pairs. The most interesting part of the phase diagram is the crossover 
regime where the two-body scattering length diverges. 

A conjectured (and, most likely, oversimplified) phase diagram for a polarized 
gas is shown in Fig. In the BEC limit the gas consists of tightly bound spin 
singlet molecules. Adding an extra up or down spin requires energy A. For \8fi\ = 
| /if — |UjJ > A the system is a homogeneous mixture of a Bose condensate and a fully 
polarized Fermi gas. One can s how that in the dilute limit this mixture is stable 
with regard to phase separation!^. 

The Bose-Fermi-gas mixture is a gapless superfluid. We can ask, along the lines 
of Sect. 13. bl whether the system is stable with respect to the formation of a non-zero 
supercurrent. We can study this question using the effective lagrangian 



Here, tp describes a gapless fermion with dispersion law e(p) and <p is the superfluid 
Goldstone mode. The low energy parameters ft and / are related to the density and 




(66) 
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Fig. 11. Schematic behavior of the gap parameter and the quasi-particle energies near the onset 
of the supercurrent state (left panel) and in the BCS limit (right panel). 

the velocity of sound. Similar to the chiral theory discussed in Sect. I3.6l the p-wave 
coupling of the fermions to the Goldstone boson is governed by the U(l) symmetry 
of the theory. 

Setting up a current v s = dip/m requires energy f 2 m 2 v 2 /2. The contribution 
from fermions can be computed using the fermion dispersion law in the presence of 
a non-zero current 



where n is the density and we have used f 2 = n/m. Son and Stephanov noticed that 
the stability of the gapless phase depends crucially on the nature of the dispersion 
law c(p). For small momenta we can write e(p) ~ e + &P 2 + Pp 4 - In the BEC limit 
a > and the minimum of the dispersion curve is at p — while in the BCS limit 
a < and the minimum is at p 0. In the latter case the density of states on the 
Fermi surface is finite and the system is unstable with respect to the formation of a 
non-zero current. The free energy functional is of exactly the same type as the one 
given in equ. (|53|l . On the other hand, if the minimum of the dispersion curve is at 
zero, then the density of states vanishes and there is no instability. As a consequence 
there is a critical point along the BEC-BCS line at which the instability will set in. 

We can also analyze the system in the BCS li mit. Th is analysis goes back to 
the work of Larkin, Ovchninikov, Fulde and Ferell 1 5^ 1 53 1 see ^h e review l^il. First 
consider homogeneous solutions to the BCS gap equation for Sfi =/= 0. In the regime 
5fi < Aq where Ao = A(5/x = 0) the gap equation has a solution with gap parameter 
A = Aq. The free energy of this solution is 



e v (p) = e(p) + v s -p-SfJ,. 



(67) 



The total free energy is 




(68) 




(69) 
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where N is the density of states on the Fermi surface. For Sfj, > A /2 there is a 
second solution with A = (25[iAq — A 2 ,) 1 / 2 , but this solution is a local maximum 
of the free energy. There are no gapless modes in the stable BCS phase, but the 
unstable BCS phase is gapless, see Fig. ITTI 

For Ao > 5 (J, > Aq/\/2 the gap equation has a non-trivial solution, but the free 
energy is higher than the free energy of the normal phase and the solution is only 
meta-stable. LOFF studied whether it is possible to find a stable solution in which 
the gap has a spatially varying phase 

A(x) = Ae 2l ^ s . (70) 

This solution exists in the LOFF window S/Ji < 6fi < 8^2 with S/ii = Ao/y/2 ~ 
0.7lAo and S/i2 — 0.754Ao- The LOFF momentum q depends on 5^. Near 6/12 we 
have qvp ~ 1.26 fi. The gap A goes to zero near 6/12 and reaches A ~ 0.25Ao at 
6 Hi. 

Clearly, the LOFF solution is of the same type as the supercurrent state. The 
U(l) of baryon number is spontaneously broken and the phase of the condensate 
has a non-zero gradient. The difference is that in the supercurrent state the current 
is much smaller than the gap, Vp (V<£>) *C A, while in the LOFF phase vp(Vtp) > A 
(and Vf(Vip) ^ A near 6^2)- In the supercurrent state the Fermi surface is mostly 
gapped but a small shell near one of the pole caps is ungapped. In the weakly 
coupled LOFF state there are gapless excitations over most of the Fermi surface 
but pairing takes place near two rings on the northern and southern hemisphere. 

In the LOFF state it is energetically favorable to take linear superpositions of 
equ. H70fl with more than one plane wave. As a consequence, the LOFF state has 
a crystalline structure and there are planes on which the order parameter has a 
node. In the supercurrent state the current is much smaller than the gap and the 
formation of nodes is not favored. This implies that there will be at least one phase 
transition (not shown in Fig. ITtj|) that separates the supercurrent state from the 
LOFF phase. 

5. Outlook 

Does QCD with three flavors more closely resemble the supercurrent state or the 
LOFF state of the cold atomic system? In the CFL phase gapless excitations first 
appear in a regime where the paired state is stable with respect to the normal 
phase or other homogeneous phases. The resulting instability can be resolved by 
the formation of a Goldstone boson current j < A. Whether or not the current is 
not only numerically but also parametrically small compared to the gap depends on 
certain details that need to be studied more carefully. In the kaon condensed phase 
the lowest mode is charged and the current is suppressed by charge neutrality. If 
the neutral mode becomes gapless, too, or if kaons are not condensed the current 
is no longer suppressed. Once the current becomes comparable to the gap it may 
be more appropriate to characterize the system as a LOFF state ^1 In particular, 
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it is p ossib le that nodes in the condensate appear and quark matter turns into a 
crystal ^3. 

Ultimately we would like to obtain observational evidence for exotic phases of 
dense matter. Polarized atomic Fermi gases have been created in the laboratory 
but so far neither the supercurrent state nor the LOFF state have been observed. 
In the case of quark matter we need to identify unique signatures of the possible 
phases that can be compared to observational evidence. Much work in this direction 
remains to be done. 
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